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ABSTRACT 

A  non-similar  boundary  layer  theory  for  air  blowing  over  a  water  layer  on  a  flat  plate 
is  formulated  and  studied  as  a  two-fluid  problem  in  which  the  position  of  the  interface  is 
unknown.  The  problem  is  considered  at  large  Reynolds  number  (based  on  i),  away  from 
the  leading  edge.  We  derive  a  simple  non-similar  analytic  solution  of  the  problem  for  which 
the  interface  height  is  proportional  to  x1/4  and  the  water  and  air  flow  satisfy  the  Blasius 
boundary  layer  equations,  with  a  linear  profile  in  the  water  and  a  Blasius  profile  in  the  air. 
Numerical  studies  of  the  initial  value  problem  suggests  that  this  asymptotic,  non-similar 
air-water  boundary  layer  solution  is  a  global  attractor  for  all  initial  conditions. 
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1  Introduction 


The  effects  of  water  layers  driven  over  solid  surfaces  by  wind  are  of  interest  in  the  performance 
of  aircraft  in  rain,  for  the  de-icing  of  airplane  wings  and  surely  in  many  other  applications. 
Since  such  problems  are  intrinsically  of  a  boundary  layer  type  and  since  the  more  interesting 
phenomenon  which  might  arise,  like  the  formation  of  waves,  film  rupture  and  the  like  are 
probably  best  framed  in  terms  of  stability,  it  is  necessary  to  derive  the  analytic  forms  that 
such  flows  will  take  on  when  instability  is  neglected.  This  derivation  is  carried  out  here. 

Previous  works  related  to  the  present  one  are  by  Yih  (1990),  who  modeled  the  de-icing 
problem,  and  by  Wang  (1992),  who  considered  the  development  of  boundary  layers  in  the 
shearing  flow  of  one  fluid  over  another.  Both  works  are  flawed  by  assuming  rather  than 
finding  the  shape  of  the  interface.  In  the  case  treated  by  Wang,  the  interface  ;s  flat  but  the 
jump  of  the  normal  stress  is  not  zero.  Yih  assumed  that  flow  in  the  air  is  a  Blasius  flow  and 
the  flow  in  the  water  is  a  simple  shear.  This  is  a  correct  form  for  the  boundary  layer,  as  we 
shall  show,  Hut  V.*.  neglects  the  variation  of  film  thickness  with  x  and  so  cannot  enforce  the 
kinematic  condition  at  the  interface  or  find  its  shape. 

In  this  work,  we  formulate  a  non-similar  boundary  layer  theory  retaining  all  terms  which 
decay  faster  than  l/£,  £  =  x(U /“Ivix)' /2,  where  is  the  kinematic  viscosity  of  air.  All  of 
the  interface  conditions  are  enforced  in  this  asymptotic  regime.  An  effect  of  the  small  but 
non-zero  vertical  velocity  component  at  the  interface  is  to  force  the  interface  to  grow  like 
x1/4  when  the  boundary  layer  in  the  air  grows  like  x1/2.  The  interface  looks  thin  on  the 
scale  of  the  boundary  layer.  Asymptotically,  at  large  x,  the  water  and  air  satisfy  the  Blasius 
boundary  layer  equations  with  a  linear  profile  in  the  water  and  the  flat  plate  profile  in  the 
air.  This  non-similar  (or  coupled  self-similar)  solution  appears  to  be  a  global  attractor  for 
all  initial  conditions. 

2  Governing  Equations 

A  water  film  of  height  y  =  h(x)  is  flowing  on  a  flat  plate  driven  by  shear  stresses  emanating 
from  an  air  stream  with  streaming  velocity  U.  We  seek  the  nature  of  the  flow  under  the 
circumstances  which  give  rise  to  Blasius  boundary  layers  in  the  flow  of  one  fluid  over  a  flat 
plate.  This  motivates  the  introduction  of  the  same  scales  that  are  used  in  the  classical  case, 
giving 
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where  u  =  /'  is  the  x  component  of  velocity,  v  =  T}f'  —  f  —  ( is  the  y  component  if  velocity, 
r  is  the  time,  p  is  the  pressure  and  /  is  the  stream  function.  The  parameters  p2  and  i/2  are 
for  air;  subscript  1  is  for  water. 

These  scales  are  introduced  into  the  Navier-Stokes  equations  which  are  written  below  in 
the  new  variables  without  approximation.  The  continuity  equation  becomes 
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The  momentum  equations  are 
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In  the  interface  equations  written  below,  [.]  designates  the  jump  (.)j  —  (.)2.  At  the  interface 
r)  =  r)*(y  =  the  velocity  is  continuous, 

[«]  =  °>  (5) 

M  =  0.  (6) 

The  shear  stress  is  also  continuous  across  t)  =  rj*, 
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The  jump  in  the  normal  stress  is  balanced  by  interfacial  tension, 
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where  G  =  ( g/U 2)  i^2v2x /U)1^  ^  g  is  the  acceleration  due  to  gravity,  S  =  T/(u2p2U),  and  T  is 
the  coefficient  of  surface  tension.  Density  and  viscosity  are  normalized  by  the  corresponding 
properties  of  the  air 
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The  kinematic  equation  of  the  free  surface  is  written  in  the  boundary  layer  coordinates 
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3  Asymptotic  solution  for  large  ( 

Here  we  exhibit  a  special  solution  of  our  coupled  air-water  system  which  will  be  shown  later 
to  be  a  global  attractor  of  non-similar  boundary  layer  solutions.  We  s.art  by  making  several 
assumptions.  First  we  assume  that  for  large  x  (large  £)  the  horizontal  velocity  component 
in  the  water  is  linear  in  y,  and  there  is  a  similarity  solution  in  the  air  with  /( 0)  =  /'( 0)  =  0, 
and  /"(0)  is  constant.  In  the  air  we  write 


and  in  the  water 


u  =  /'(*)»  *  = 

c(x). 
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u  =  c(x)y,  Q  =  /  udy  =  ~r^h2  =  constant. 
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Secondly  we  assume  that  the  interface  position  77*  — ►  0  as  £  — ►  00.  This  assumption  is 
equivalent  to  assuming  that  as  £  — ♦  00  the  liquid  layer  is  a  vanishingly  small  fraction  of  the 
boundary  layer  in  the  air. 

We  are  going  to  show  that  the  continuity  of  the  shear  stress  across  the  interface  77*  implies 
that  c  =  Ax-1/2,  with  A  a  constant,  so  that  u  =  At]  can  be  expressed  in  terms  of  77  alone. 
The  interface  is  on  y  =  h(x),t]  =  77*  =  k-[ 75-.  The  shear  stress  is  continuous 
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Using  (13),  we  may  write 


c(x)  = 


du 

12-K-  at  77 
dy 

1  ]Ty- 

7*  =  h(x) 

M/"(  0) 


fiix 


1/2 


2  Q 
A  ' 


x'!r 

(14) 

k  ...  kh  . 

^x'/*f  V/2'' 

(15) 

(/x1^2  — ►  0,  so  that 

— ►  /"( 0)  which 

A 

i1/2 

(16) 

BV/J, 

(17) 

3 


h  =  Bx'/'. 
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Using  (18),  we  find  that 


Now 
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For  large  x,i )*  — >  0,  so  that  f(r)*)  — ►  /'( 0)  =  0.  Actually  we  know  that  u  tends  to  zero 
because  the  mass  flux  is  constant  but  the  flow  area  (18)  keeps  increasing.  We  may  also  show 
that  v(h(x)yx)  tends  to  zero  for  large  x.  We  can  do  this  in  different  ways.  Suppose  we 
multiply  the  kinematic  condition  (11)  by  £u  and  let  £  — ►  oo. 
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From  (19)  we  know  that,  since  £  =  kx1?2,  as  £  — >  oo, 
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Thus  as  £  — ►  oo, 

dr)*  k3>2B  n 
*  d£  ~  £*/2  ->0, 

Since  u  — >  0  as  £  — *  oo,  so  also  v  — ♦  0.  So  as  £  — >  oo,  both  u(h(x),  i)->0  and  v(h(x),  x)  —+  0; 
that  is  f(rj*)  — *  0  and  /(??*)  — *  0.  The  interface  conditions  then  become  the  same  as  the 
boundary  conditions  for  a  Blasius  boundary  layer  in  the  absence  of  the  water  layer  but 
displaced,  and  so  the  velocity  in  the  air  is  given  by  the  Blasius  solution  for  the  boundary 
layer  over  a  flat  plate.  The  solution  in  the  water  f  =  At)  also  satisfies  the  Blasius  equation 
/  f"  +  vi  f"  =  0,  0  <  r)  <  r)*  — >0  and  implies  that  the  jump  conditions  [/]  =  [/']  =  0  may 
be  replaced  with  the  Blasius  boundary  conditions  /( 0)  =  /'( 0)  =  0  applied  upon  a  solid 
surface. 

This  shows  that  our  original  assumptions  are  self-consistent.  The  reader  may  verify 
that  the  solution  (13)  with  c(x)  given  by  (16)  is  self-similar  (satisfies  (24)  without  the  the 
bracketed  term).  The  main  good  luck  that  we  have  is  that  r)*  — ►  0  for  large  x  even  though 
h(x)  — >  oo 


4  Non-similar  boundary  layers 

Since  the  asymptotic  solution  for  large  £  is  similar,  but  not  in  the  usual  sense,  we  are 
motivated  to  see  if  this  solution  can  be  embedded  in  a  large  class  of  non -similar  solutions  for 
large  £  (large  Reynolds  number  Ux/ 2u2)  and  to  investigate  the  possibility  that  the  asymptotic 
solution  enjoys  a  special  status  as  a  global  attractor  for  all  solutions  of  this  non-similar  family. 

Self  similar  boundary  layers  depend  only  on  tj  and  not  on  £.  The  £  derivatives  of  u,  v  and 
Tj*  may  be  small,  of  the  order  of  £-n,n  >  1,  when  £  is  large.  We  retain  those  terms  in  each 
equation  which  are  0(1),  dropping  all  terms  of  0(£-n),  For  this  computation  we  assumed 
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that  the  first  derivative  of  u  with  respect  to  £  scale  with  <f  1 ,  the  second  derivative  with  £  2 
and  =  0(r3/2)-  These  scalings  can  be  verified  a  posteriori.  We  find  that 


/r + - 


c,r^L_r<tL)  =  o 

s(;  ^  j  df) 


(24) 


and  at  r?  =  77*: 

|j  =  0,  [/1  =  o,  f/j  =  o,  (></"]  =  o,  W  =  Wot”.  ^  =  (**) 

Equation  (24)  can  be  found  in  Schlichting  (1987).  We  may  remark  that  the  contributions 
of  the  viscous  terms  in  the  normal  stress  balance  are  £~2  times  the  terms  retained  and  the 
surface  tension  terms  are  £-3  times  those  retained. 
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5  Solutions  of  initial  value  problems  for  the  non-similar, 
two-fluid  boundary  layer  equation 

5.1  Evolution  of  a  profile  which  is  initially  parabolic  in  segments. 

Since  the  non-similar  solutions  depend  on  £,  we  are  obliged  to  consider  the  evolution  of  flows 
prescribed  at  some  initial  position  £  =  £0.  Fortunately  our  simulations  of  those  initial  value 
problems  indicate  a  rapid  decay  from  initial  values  to  the  asymptotic  solution  described  in 
§3.  This  situation  is  not  unfamiliar.  In  the  classical  theory  of  boundary  layers,  the  Blasius 
solution,  the  similarity  solution,  is  a  global  attractor  for  all  initial  value  problems  which  are 
not  similar  (see  Serrin  (1967)).  In  our  case  the  attractor  cannot  be  self-similar,  but  the  non¬ 
similar  solution  of  §3  arises  asymptotically,  for  large  £  from  all  initial  conditions  explored. 
The  first  initial  condition  we  used  to  solve  (24)  and  (25)  is  plotted  below. 


(a)  (b) 

Figure  1.  Parabolic  initial  condition  assumed  for  the  horizontal  velocity  component  /'  with 
Tj*  =  0.15  and  £0  =  50.  (a)  Initial  condition,  (b)  Detailed  plot  of  initial  condition  at  the 
interface. 

For  this  initial  condition  77*  =  0.15  and  £<,  =  50.  These  values  were  chosen  to  be  representative 
of  conditions  discussed  by  Hastings  and  Manuel  (1985),  in  which  they  describe  the  results 
of  their  wind  tunnel  measurements  of  a  wing  in  simulated  rain.  In  the  above  profile,  for 
0  <  r;  <  i]* ,  f  is  parabolic  with  f"  being  positive.  For  t;*  <  rj  <  3.5,  /'  is  parabolic  with  /"' 
being  negative,  and  for  J]  >  3.5,  f  —  1.0. 

We  solved  (24)  and  (25)  subject  to  the  aformentioned  initial  conditions  for  >  <f0  using 
a  finite  difference  scheme  found  in  Schlichting  ( 1987).  This  scheme  is  iterative  using  second- 
order  differencing  for  the  tj  derivatives,  and  a  first-order,  implicit  Euler  differencing  for  the  £ 
derivatives.  We  assumed  an  interface  r/*  =  a£u,  where  a  =  (to)  /  aiu^  found  the  n(£) 
given  in  figure  2. 
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Figure  2.  Plot  of  «  for  the  non-similar  boundary  layer  with  initial  conditions  from  figure  1. 


The  exponent  function  n(£)  decays  rapidly  to  -1/2,  corresponding  to  a  limiting  power 
law  increase  h(x)  <x  x1^4  of  the  interface  height  as  given  by  the  asymptotic  solution  shown 
in  figure  1.  For  £  not  much  larger  than  £0  the  horizontal  velocity  component  /'  becomes 
linear  in  ?/,  /'  =  c(£)r],  in  the  water  layer  and  remains  so  for  all  greater  £.  The  shear  stress 
in  the  water,  is  directly  related  to  the  shear  stress  in  the  air  at  the  interface  through 
the  tangential  stress  condition  found  in  (25).  From  our  calculations  we  find  that  the  shear 
stress  in  the  air  at  the  interface,  /",  asymptotically  tends  to  the  shear  stress  in  the  Blasius 
boundary  layer  of  a  single  fluid  over  a  flat  plate  at  the  plate  surface,  as  shown  in  figure  3. 


/'<*  >0*)  k 


(a)  (b) 

Figure  3  (a)  Detailed  plot  of  boundary  layer  solution  in  the  water  layer  at  £  =  50.0675  for 
the  initial  conditions  of  figure  1.  (b)  Shear  stress,  /",  in  the  air  at  the  interface,  77*  =  a£“'^, 
for  £  > 

Notice  that  the  exponent  function  n(£)  tends  to  its  asymptotic  value  extremely  rapid,  while 
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the  shear  stress  at  the  interlace  in  the  air  /"  has  a  much  slower  decay.  The  reason  for  this 
difference  is  n(f)  becomes  nearly  -1/2  as  soon  as  the  horizontal  velocity  component  in  the 
water  becomes  linear,  which  occurs  for  if  not  very  much  greater  than  f„.  Figure  3  shows  that 
the  How  has  not  vet  reached  its  as'c  ptotic  form,  since  the  shear  stress  at  the  interface  in 
the  air  is  appreciably  different  f:  i  its  asymptotic  value.  The  flow  then  slowly  tends  to  its 
asymptotic  state.  The  flow  b  on-similar  for  finite  f  because  the  boundary  layer  height  and 
the  interface  do  not  seal*3  the  same.  Asymptotically  one  may  say  the  flow  is  similar  because 
i/*  goes  to  zero  and  the  flow  in  the  air  is  the  self-similar  Blasius  profile. 

5.2  Evolution  of  another  profile  which  is  initially  parabolic  in 
segments. 

In  the  next  example  we  used  the  same  initial  /'( rj )  for  T]  >  r?*  as  in  figure  1(a),  but  the  profile 
shown  in  figure  1(b)  was  replaced  with  the  one  shown  in  figure  4. 


Figure  4.  Initial  profile  in  the  water  at  f  =  The  profile  in  the  air  is  shown  in  figure  1(a). 

In  figure  5  we  show  the  evolution  of  the  exponent  n(f)  for  f  >  f 0  for  the  aformentioned 
initial  condition.  Obviously  n(f)  — >  —1/2  for  large  f  and  77*  =  af~1/2  for  large  f.  The 
solution  given  in  §3  is  attained  asymptotically. 
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5.3  Evolution  of  a  profile  which  is  initially  in  similarity  form  but 
with  tj*({0)  =  0.15  >0. 

The  initial  condition  shown  in  figure  6  is  generated  from  (24)  and  (25)  when  the  £  derivatives 
are  set  to  zero,  so  that 

//" + (26) 

and  at  rj  =  r}m: 

^  =  0,  [/']  =  0,  (/]  =  0,  \nf")  =  0,  \p\  =  [pW.  (27) 

This  system  is  solved  using  a  iterative  finite  difference  scheme  using  second  order  differencing 
and  gives  rise  to  the  profiles  shown  in  figure  6. 


n 


(a)  (b) 

Figure  6.  Solution  /'(tj)  of  (26)  and  (27).  (a)  rj  >  0.15.  (b)  //  <  0.15.  The  Blasius  solution 
of  (24)  for  /( 0.15)  =  // (0.15)  =  0  is  also  shown  in  (b). 
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The  profiles  in  figure  6  were  used  for  initial  values  in  the  non-similar  equations  (24)  and 
(25)  and  they  evolved  to  the  coupled  self-similar  solutions  described  in  §3  in  which  t?*(£)  — ♦  0. 
The  difference  between  the  solution  in  the  air  and  the  classical  Blasius  flat  plate  solution 
also  tends  to  zero.  A  graph  of  the  exponent  function  for  this  example  is  shown  figure  7. 


6  Behavior  of  £  derivatives  of  the  non-similar  solution 

In  deriving  the  system  which  lead  to  the  above  solutions,  we  assumed  that  first  and  second 
derivatives  of  u  and  rj*  with  respect  to  £  were  inversely  proportional  to  £“",n  >  1.  ^From 
our  solution  we  find  that  ^  scales  like  £~3'2  and  scales  like  £~5/2,  in  agreement  with 
our  scalings.  Since  £  scales  with  a:1/2,  the  interface  position  scales  with  x1/4. 

In  order  to  examine  the  £  derivatives  of  u,  we  define  two  functions, 

F,M)  =  f§p  (28) 

(29) 

Depending  upon  the  initial  conditions  chosen  at  £0,  Fi(rj,£0)  and  F2(  77,  £0)  may  be  extremely 
large.  However,  for  larger  values  of  £  the  horizontal  velocity  component  in  the  water  becomes 
linear  in  r/,  n(£)  tends  to  -1/2  and  Fi(t;,£)  <  1,  F2(77,£)  <  1  (see  figure  8). 
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Figure  8.  Plot  of  Fi(rj,()  and  versus  rj  for  £  =  50.01  for  the  flow  plotted  in  figure  6. 

|Fj|  and  IF2I  decrease,  as  shown  in  figure  9,  although  the  decrease  in  |Fij  is  very  small. 
Since  as  £  — ►  00,  tj*  — ►  0  and  the  flow  in  the  air  goes  to  the  Blasius  boundary,  both  F\  — *■  0 
and  F2  — >  0  as  £  — ►  00. 


Figure  9.  Plot  of  F\{tf ,£)  and  F2(t7£)  versus  77  for  (  =  129.1  for  the  flow  plotted  in  figure  6. 


Then  since  F\  and  F2  are  both  0(1),  the  assumptions  that  were  made  in  deriving  (24)  and 
(25)  are  shown  to  hold. 


7  Concluding  remarks 

The  analysis  given  in  this  paper  can  be  extended  to  two-fluid  boundary  layer  problems  with 
other  free  streams,  say  U  =  UooXn.  It  is  also  probable  that  the  solution  given  in  §3  and  the 
other  limiting  solutions  to  which  we  have  just  alluded  are  unique  large  x  limits  of  steady 
coupled  air-water  solutions  of  the  Navier-Stokes  equations  with  different  initial  profiles  at 

X  —  X  Q  . 
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